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We apply local orbital basis density functional theory (using SIESTA) coupled with a mapping to 
the Anderson impurity model to estimate the Coulomb assisted or correlated hybridization between 
transition metal d-orbitals and ligand sp-orbitals for a number of molecular complexes. We find 
remarkably high values which can have several physical implications including: (i) renormalization 
of effective single band or multiband Hubbard model parameters for the cuprates and, potentially, 
elemental iron, and (ii) spin polarizing molecular transistors. 

PACS numbers: 31.10,+z, 71.15.-m, 71.27.+a 
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The properties of transition metal compounds are of- 
ten dominated by the 3 c? orbitals because their localized 
character causes a strong Coulomb repulsion between the 
electrons. The theoretical progress in the field has been 
impeded by the extreme difficulties of dealing with even 
the simplest model Hamiltonian appropriate for these 
materials, the Hubbard model consisting of a direct 
electron hopping between orbitals i and j with amplitude 
Uj and of the Coulomb repulsion U of electrons in the 
same orbital. U is a matrix element of the Coulomb po- 
tential, U = (ii\V(r — r')\ii). Matrix elements involving 
different lattice sites i and j are generally smaller than 
U. One of them is the so-called correlated hybridization 
Xij = (ii\V(r — that describes a density depen- 

dent hopping. Estimates for (0.5 eV in transition 
metals p|, 0.8 cV in cuprates @|) show that the matrix 
elements of the correlated hybridization are comparable 
to or even larger than the corresponding amplitudes of 
direct hopping in real systems. Up to now, effects caused 
by correlated hybridization have not yet attracted much 
attention: Hirsch suggested a new mechanism for super- 
conductivity 0, some exactly solvable cases were dis- 
cussed 0, finite clusters were studied using exact diag- 
onalization |j| , metallic ferromagnetism was investigated 
, and correlated hybridization was studied via dynam- 
ical mean field theory in the Falikov-Kimball model Q. 

In many compounds, cuprates and manganites are 
well-known examples, the transition metal atoms are sur- 
rounded by oxygens or other elements with p-orbitals so 
that the 3d orbitals are only effectively coupled with each 
other due to oxygen orbitals. Realistic multi-band mod- 
els including oxygen degrees of freedom, as proposed for 
the cuprates jg, can not easily be mapped onto effec- 
tive single-band models (see j!| and references therein). 
The effective hopping between transition metal sites must 
then be mediated by the hybridization with p-orbitals of 
the oxygens. Despite the somewhat robust interest in 
correlated hybridization between like orbitals, there have 
been no ab initio studies of the corresponding correlated 
hybridization |ldl | between the d orbitals and surround- 
ing p orbitals. This may be of considerable interest in 
molecular transistors based upon transition metal com- 



plexes for which the surrounding ligand atoms provide 
the linkage to the leads [lll| . 

In this paper, we estimate matrix elements of the corre- 
lated hybridization for several transition metals in differ- 
ent chemical environments from density functional theory 
(DFT) calculations. To our knowledge, this is the first 
systematic ab initio study of correlated hybridization ma- 
trix elements, and we find remarkably high values, al- 
though similar spin dependent hybridization phenomena 
have been noted in the context of molecular transistors 
[T^ . We demonstrate that these additional hybridization 
matrix elements could significantly change the parame- 
ters in effective single band models for transition metal 
oxides. Furthermore, we demonstrate that the correlated 
hybridization may possibly provide a means to significant 
spin polarization of currents through transition metal 
based molecular transistors in modest magnetic fields. 

We have carried out spin-polarized electronic structure 
calculations using the fully ab initio DFT code SIESTA 
|l3l |. It uses Troullier-Martins norm-conserving pseudo 
potentials [l4| in the Kleinman-Bylander form where 
we included nonlinear partial-core corrections for the 
transition metal atoms to take into account exchange 
and corelation effects in the core region [l^. We used 
the generalized gradient approximation (GGA) for the 
exchange-correlation energy functional in the version of 
Ref. hj. SIESTA uses a basis set of atomic orbitals 
where the method by Sankey and Niklewski [l^ is em- 
ployed. We used a double-^ basis set, and included po- 
larization orbitals for the transition metal atoms. To de- 
termine the minimum energy configuration within DFT 
all complexes were allowed to relay till the force on each 
atom was less than 0.03 eV/A. 

The results of the SIESTA calculations are used to 
determine effective hybridization matrix elements. For 
that purpose we map the final DFT Hamiltonian 7Ydft 
obtained from SIESTA onto an effective two-band model, 

w = ^4,4*4* + '52% ur pLpc U , (i) 

+ (4x?w + h. 



TABLE I: Maximal difference between the spin directions 
of the effective hybridization matrix elements Vi aa for sev- 
eral transition metal complexes in the range ±10 eV about 
the chemical potential where the overlapp between the ligand 
states of the two spin directions is at least 0.95. A max is de- 
fined in Eq. @ an d the last column gives the orbital symme- 
try of the maximally different effective hybridization matrix 
elements referenced to tetrahedral or octahedral symmetry. 



Compound 


tmx(\V ia i - V ial \) 




Orbital 


[Fe (SH) 4 ] 2 ~, perfect structure 


0.361 eV 


0.791 


h 9 


[Fe (SH) 4 ] 2_ , relaxed structure 


0.481 eV 


0.354 


h 3 


[Ou (H 2 0)J 2+ 


0.871 eV 


0.875 


e s 


[Mn (H 2 0) 6 f + 


2.416 eV 


2.000 


h 3 


[Ni (H 2 0) 6 ] 2+ 


1.584 eV 


1.482 




Co(3,5-DBSQ) 2 (phen) 


1.021 eV 


4.700 


h s 


Co(3,5-DBSQ)(3,5-DBCat)(phen) 


0.271 eV 


0.542 


h a 


Mn(N0 2 -phen) (3,6-DBSQ) 2 


0.743 eV 


0.664 


h 3 


Mn(II) (Bupy) 2 (3,6-DBCat) 2 


0.479 eV 


0.555 


U g 


Mn(III)(Bupy) 2 (3,6-DBCat) 2 


0.987 eV 


0.717 


e 9 


Mn(IV) (Bupy), (3,6-DBCat) 2 


0.894 eV 


1.132 





where the 3d orbitals of the transition metal atoms are 
separated from the rest. The d irJ (p arT ) are fermionic cre- 
ation operators of 3d (ligand) electrons with spin a at 
orbital i (a). To determine the onc-particlc energies, ef a 
and S^ a , and the hybridization matrix elements, Vi aa , 
we apply unitary transformations to TLbft- Because the 
atomic basis states used by SIESTA are not orthogonal 
the transformation of TLdft is performed in three steps: 
(i) The block of the 3d orbitals (single £ only) is diagonal- 
ized. (ii) The 3d contributions to the ligand basis states 
are removed, (iii) The ligand block is diagonalized. Af- 
ter theses steps we obtain an effective Hamiltonian TL of 
the desired form where TL has the same eigenvalues 
as 7Ydft because TL is derived by unitary transformation 
from TLdft- 

Because we performed spin-dependent DFT calcula- 
tions TLdft consists of two completely separated spin sec- 
tors. Therefore, the parameters of TL depend on the spin 
direction, and we find a remarkable spin dependence of 
the hybridization matrix elements if the transition metal 
atoms are in high-spin states (see Tab.[Q|. 

To discuss the spin dependence of the hybridization 
matrix elements in detail we calculate the hybridization 
broadening function, T ia (uj) = J2 a \V laa \ 2 5(uj - e p acr ), for 
all impurity states i. As one can see from Fig. ^ the 
hybridization matrix elements of the minority spin direc- 
tion are always larger than the values of the majority 
spin direction. 

A detailed analysis shows that this spin dependence 
can only derive from correlated hybridization. To illus- 
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1: (Color online) Spin dependent hybridization func- 

1 2- 



tions for [Fe (SH) 4 J , perfect structure, where the line spec- 
tra have been broadened with Gaussian functions of width 
0.5 eV . The results for the majority [minority] spin direction 
is shown in panel (a) [panel (b)]. 



trate the point, let us now consider a Hamiltonian, 

« = E £ i4^+5>?>Uw (2) 



d\ a Paa + H.C. 



with correlated hybridization matrix elements Ai a . TL 
can be easily transformed to the form of Eq. Q if a 
factorization approximation is used where the correlated 
hybridization leads to spin dependent one-particle ener- 
gies if a and hybridization matrix elements Vi aa . This 
consideration can also be used to determine the matrix 
elements A; Q of the correlated hybridization, one obtains 



A, 



Vi 



Vial 



(3) 



Because of the expectation values of the occupation num- 
ber operators are restricted, < (d icr di„) < 1, the dif- 
ference between the both spin directions of the effective 
hybridization matrix elements (as listed in Tab. [IJ has 
to be interpreted as a lower bound for the full matrix 
element of the correlated hybridization. Thus, we con- 
clude from Tab. [I] that the correlated hybridization ma- 
trix elements Ai a are comparable in magnitude to the 
regular tight-binding hopping amplitudes Vi a . Although 
the commonly used approximations of DFT such as local 
density approximation or CGA are known to fail in ac- 
curately describing strongly correlated transition metal 
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systems (notably the cupratcs and metallic plutonium), 
there is a notable tradition of successfully using DFT to 
provide estimates of parameters for many-body Hamilto- 
nians, as in calculations of the hybridization for success- 
fully estimating Kondo scales in a series of cerium heavy 
fermion systems |Tflj . This is the spirit of our approach; 
while we do not include the feedback modification of the 
density to the DFT, that has successfully been done only 
in one instance pcj . 

In the following we want to show that the large val- 
ues obtained for the correlated hybridization matrix el- 
ements could indeed lead to intriguing physical effects. 
For this purpose we consider in the following with the 
cuprates and a molecular device two physical systems 
where the local coordination environment of the transi- 
tion metal atoms is precisely the same as in the clusters 
studied above. At first, we consider a three-band model 
as proposed for the cuprates Q where we use standard 
parameters [2l|, Cu 3d O 2p hybridization t p d = 1.3 eV, 
O 2p O 2p hybridization t pp = 0.65 eV, Cu on-site repul- 
sion Ud = 8.8 eV, and charge transfer energy A = 3.5 eV. 
Furthermore, we add a correlated hybridization between 
Cu 3d and O 2p orbitals to the Hamiltonian with matrix 
elements A p d = 0.5 cV. In accordance to our DFT esti- 
mates, A p d has the same sign and phase factors as t p d- 
The three-band model can be reduced j2^| to a S = | 
Heisenberg model on the square lattice of Cu sites where 
the coupling strength in fourth order perturbation theory 

is given by J = + <^+ A ^ 2 Thus, the matrix el- 
ements A p d = 0.5 eV of the correlated hybridization lead 
to an increase of the Heisenberg exchange from 0.24 eV 
to 0.34 eV. 

The correlated hybridization also affects the mapping 
of the three-band model onto an effective single-band 
model by means of a cell perturbation theory 0, l23j . 
Despite an increase of the magnitudes of the effective 
hopping matrix elements by about 20%, the on-site Hub- 
bard U is significantly lowered from 2.48 eV to 1.94 eV. 
This reduction in U may have some relevance in under- 
standing the need to use smaller Hubbard interactions in 
models of, e.g., metallic Fe than one obtains from esti- 
mates using constrained occupancy DFT [24| : If there is 
significant correlated hybridization of d levels with s,p 
levels on a neighboring site then the same mechanism of 
U reduction upon folding down to the d-band only model 
from the multi-band model will be operable. 

Perhaps the most interesting application of our re- 
sults is to tunable molecular transistors 0, |2f| |2(| , with 
the strongest response possible not in the Kondo limit 
but rather the mixed valent limit. In the presence of a 
relatively modest field, the correlated hybridization can 
impart a potentially significant spin polarization of the 
current through the device without requiring spin polar- 
ization of the leads as has been considered both exper- 
imentally [23 and theoretically [2^]. To illustrate the 
idea, consider the schematic device of Fig. [21 m which a 
transition metal based molecule is attached to ordinary 
metallic leads, and assume that a spin 1/2 large U model 
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FIG. 2: Panel (a) shows the schematic design of the proposed 
spin filter device with the definitions of the lead-to-ligand and 
ligand-to-transition metal site hybridization matrix elements 
W and V a - The energies of the involved electronic states are 
sketched in panel (b). 



is applicable, as relevant for the low spin Co complexes 
considered in Ref. |lll | (electron hopping) or a Cu com- 
plex (hole hopping). The bare lead-to-transition metal 



site hybridization Vl<i goes as V[ d w si-e f wnere W 
is the lead-to-ligand hybridization matrix element, £l is 
the ligand LUMO (HOMO) for electrons (holes) to tran- 
sition from lead to ligand, and V a = Vq(1 — Amu) is 
the spin dependent hybridization discussed earlier, with 
A a dimensionlcss ratio of the assisted hopping matrix 
element to the direct one Vq, and m the induced polar- 
ization m = rid-\ — nd{ of the transition metal ion in an 
applied field H . Note that estimates, 



A n 



= 2- 



max I Vi 



ia J 



Via 



Vial + Vial 



(4) 



for the ratio A can be found in Tab. [I] Importantly, we 
note that A can be as large as 3-5. 

We next assume that the Friedel sum rule (2j| can be 
applied so that the number of electrons (holes) rid on the 
transition metal site is given by rid = + n di = T ^ 1 
and the polarization m by m = — (6f — 8^) where 8^ 
is the Fermi energy phase shift for an electron of spin 
a to scatter off the transition metal site. We further 
take a renormalized resonant level model for the phase 
shift, with resonance position e a = eq — fiHa, \i the 
magnetic moment of the d-levels, and resonance width 
T a = To(l — 2Ama) where only contributions linear in 
the magnetic field H are taken into account. Note that 
in the Kondo regime the zero field resonance parameters 

satisfy \J Eq + f § w ksTn, with Tk the Kondo temper- 
ature of the transistor. In this model, the phase shift 
is 8 u = j — tan -1 (f 2 -)- Assuming zero bias, iden- 
tical leads, and low temperature, the conductance Q a 
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within our mean-field treatment of the correlated hy- 
bridization in spin channel a is approximately given by 
HI, G a « ^sin 2 <5 a . 

Using the definitions and results of the previous para- 
graph we thus find 

(g T -gj) 2 S in(2£ ) 

Q(H = 0) ~ i + M sin(2 5 ) f 1 ' 

with Sq the zero field phase shift. 

The A dependent enhancement factor follows from the 
impurity polarization m oc — . This effect is overes- 
timated within our mean field approximation, although 
it has a clear physical origin: the up-spin resonance nar- 
rows significantly in applied field as the down spin reso- 
nance broadens, allowing a feedback effect that enhances 
the tendency towards magnetization saturation. The ef- 
fect is most pronounced in the mixed valent regime for 
rid ~ 1-5 where a critical A value of 7r/2 yields divergent 
m within the approximation. Indeed, a significant en- 
hancement of the zero field local susceptibility was found 
in a nonperturbative treatment of this model for Uj = 
in the mixed valent regime [Tof. It might be possible 
to tune the resonance into the highly polarizing regime 
through a combination of gate control (to tune e^) and 
contact chemistry (to tune the lead-transition metal ion 
hybridization). The enhancement effect can also be sig- 



nificant in the Kondo regime, for which 1 — nd\ < 0.3 [3ll | 
and | sin(2#o)| < 0.81. Examining the A values in our ta- 
ble, we obtain critical |1 — nj\ values of 0.11 and 0.29 for 
the Co(II)(3,5-DBSQ) 2 (phen) and [Mn(H 2 0) 6 ] 2+ com- 
plexes, safely within the Kondo regime. 

To conclude, we have presented here the first system- 
atic ab initio study of correlated hybridization matrix 
elements for several transition metal complexes. Based 
on DFT calculations we have found correlated hybridiza- 
tion matrix elements comparable in magnitude to the 
regular tight-binding hybridization amplitudes which can 
lead to significant changes in the parameters of effec- 
tive models for transition metal compounds. Finally, 
we have sketched how the correlated hybridization can 
be employed to design spin sensitive devices, polarizing 
the spin current dramatically with potentially modest ap- 
plied fields and paramagnetic leads. 
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